Solitons in quadratic nonlinear photonic crystals 
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We study solitons in one-dimensional quadratic nonlinear photonic crystals with modulation of 
both the linear and nonlinear susceptibilities. We derive averaged equations that include induced 
cubic nonlinearities and numerically find previously unknown soliton families. The inclusion of the 
induced cubic terms enables us to show that solitons still exist even when the effective quadratic 
nonlineanty vanishes and conventional theory predicts that there can be no soliton. We demonstrate 
that both bright and dark forms of these solitons are stable under propagation. 
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The physics and applications of photonic band-gap 
(PBG) materials, or photonic crystals, have been active 
topics of research for more than a decade. The theory 
of linear photonic crystals is now well understood, and 
many of their fundamental properties and technical ap- 
plications have been characterized The next impor- 
tant step in the application of photonic crystals is to cre- 
ate tunable PBGs. Tunability is possible in linear pho- 
tonic crystals through, e^., the temperature dependence 
of the refractive index g or the electro-optic effect [||. 
Ultrafast dynamical tunability of the PBG can be ac- 
complished using nonlinearity, as was first demonstrated 
with a constant Kerr nonlinearity 

Here we consider quadratic nonlinear photonic crystals 
(QNPCs) that have a linear grating (periodic dielectric 
constant) and/or a nonlinear grating (periodic second- 
order or x^^^ susceptibility). QNPCs are of interest for 
all-optical components due to the fast and strong nonlin- 
earity they can provide through the parametric cascading 
effect 1^ . The efficiency of the cascading process depends 
critically on the phase mismatch between the fundamen- 
tal and second-harmonic (SH) waves, but two powerful 
methods exist that use exactly a periodic photonic crystal 
structure to control the mismatch [||-^ . In one method, 
a QNPC with a linear Bragg grating is used to bend the 
dispersion curve near the PBG [fIjlO). However, the short 
period, which is of the order of the optical wavelength, 
can be inconvenient. The second scheme, quasi-phase- 
matching (QPM), controls the phase mismatch us- 
ing a nonlinear grating with a period equal to the com- 
paratively long beat length (typically of the order of mi- 
crons). QPM is also possible with linear gratings, but 
this is much less effective [^]. 

One of the spectacular manifestations of nonlinearity 
is the soliton, a self-localized entity that can propagate 
unchanged over long distances. Homogeneous x*"^"* 
terials support solitons in all dimensions and gap 
solitons exist in QNPCs with a linear Bragg grating p2[ . 
In this Letter we focus on the open fundamental problem 
of whether solitons exist in ID QNPCs with both a linear 
and a nonlinear QPM grating. Such a simultaneous linear 
grating is difficult to avoid when, for example, creating 
nonlinear QPM gratings in GaAs/AlAs semiconductors 
through quantum- well disordering jl^ ]. 



Solitons exist in ID QNPCs with a nonlinear QPM 
grating , but a simultaneous linear grating can reduce 
the effective x*-^"* nonlinearity [0-1). Thus the global ex- 
istence of solitons in such QNPCs is nontrivial. Building 
on previous findings that nonlinear QPM gratings induce 
cubic nonlinearities (m) , we find soliton solutions that are 
stable under propagation. The induced cubic terms pro- 
vide an elegant means of showing that, against intuition, 
the QNPC supports stable bright and dark solitons even 
when there is no effective x^^-* nonlinearity. This is anal- 
ogous to the existence of solitons in dispersion-managed 
fibres with no average dispersion [|l5|. 




FIG. 1. Geometry of the ID QNPC in the form of a x^^^ 
slab waveguide. Black and white domains indicate regions 
with different dielectric constant and x''^'' coefficient. 

We consider the interaction of a cw beam (carrier 
frequency uj) with its SH, propagating in a lossless 
ID QNPC under conditions for type I second-harmonic 
generation (SHG), as sketched in Fig. |l|. We assume 
that the modulation of the refractive index is weak 
{Anj(z)/fij<^l, where nj{z)=nj + Anj{z) and j refers 
to the frequency juj), and we consider only gratings for 
forward QPM. The grating period is then much longer 
than the optical period, in which case Bragg reflections 
can be neglected. The evolution of the slowly varying 
beam envelopes is then described by pjlq] 
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where w=w{x, z) and v—v{x, z) are the envelope func- 
tions of the fundamental and SH, respectively. The trans- 
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verse and propagation coordinates x and z are in units 
of the input beam width xq and the diffraction length 
Ld=kixl, respectively. The parameter f3=AkLd is pro- 
portional to the mismatch Ak=k2-2ki, kj=jujnj / c be- 
ing the average wavenumber. Thus [3 is positive for 
normal dispersion and negative for anomalous disper- 
sion. The normalized refractive index grating is given by 
aj{z)=LdLoAnj{z)/c and the normalized nonlinear grat- 
ing by x{z)^LdUjdcs{z) / {fiicV where d^s^x^'^V'^ is given 
in MKS units. The model (P describes both temporal 
and spatial solitons fi^. 
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FIG. 2. Normalized linear and quadratic nonlinear grat- 
ings, aj{z) and x{z)j with period 2Lo=2Tv/\Hi\. 

The aim is now to average Eqs. and derive accurate 
equations for the average field. To do so we focus on first 
order QPM using the conventional square gratings with 
50% duty cycle, shown in Fig. ||. We expand the grating 
functions in Fourier series: 
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where gn—2s/{iTm) for n odd and gn=0 for n even, with 
s = sign(K). The gratings drive the system, which means 
that we may expand the envelope functions in Fourier 
series also: 



w = 2_^Wn{z,x)e'"'''', V 



^^;„(^,x)e^("--'^)^ (3) 



assuming that the coefficients Wn{z,x) and Vn{z,x) vary 
slowly in z compared to exp(iKz). The residual mismatch 
f3=f3 — K is ideally zero. 

Three physical length scales are in play: the diffrac- 
tion length Ld, the coherence length Lc, and the grating 
domain length Lq. In normalized units Ld—1, Lc—7r/\f3\, 
and Lo=tt/\k\. We assume a typical QPM grating with 
a domain length that is much shorter than the diffrac- 
tion length, io^Cl. Furthermore, the grating is of good 
quality, with the domain length being close to the coher- 
ence length, Lq—Lc^ so the residual mismatch is small, 
|/3|<C|/v|. In this case and we can use perturbation 

theory with the small parameter e=l/jK\ ^ 1. 

Following the approach of Ref. |jl4], we insert the 
Fourier expansions and (^) into the dynamical equa- 
tions and assume the harmonics Wn^o and Vn^o to be of 
order e. To lowest order (e^), this gives the harmonics 



Wn^o = {aignWo + dgn-iWQVo)/{nK), 
Vn^o = (2a2g„wo + dgn+iwl)/{nK). 



(4) 



Using these solutions, we obtain to first order e the aver- 
aged equations for the DC components wq and vq'. 

.dwo Id^Wo , , ,2 I |2n n 

2~dx^ P^^o^o + 7(l^^o| " \wo\ )wo = 0, 
~dz 4:~dx^ + pvo+ p Wq+ 2j\wo\ vo = 0. (5) 

These equations also describe mth order QPM (where 
f3=f3 — rriK is ideally zero) and any other type of peri- 
odic grating, the parameters p and 7 being simply given 
as sums over the Fourier coefficients of the grating . 
Incorporating time or the spatial y coordinate is also 
straightforward. For the square grating (|2|), p and 7 can 
be explicitly calculated to 
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From Eqs. follows the important result that cubic 
nonlinearities are induced in QNPCs by nonlinear QPM 
gratings. This cubic nonlinearity has the form of self- 
and cross-phase modulation (SPM and XPM), and is a 
result of non-phase-matched coupling between the wave 
at the main spatial frequency k and its higher harmon- 
ics. It is thus of a fundamentally different nature than 
the material Kerr nonlinearity, which is reflected in the 
fact that the SPM term is absent for the SH. 

The averaged model (||) is identical to the known 
model for nonlinear QPM gratings with no DC compo- 
nent {aj=do=0), for which simulations confirmed that 
bright solitons had properties that were not predicted by 
the conventional model with only quadratic terms, but 
were accurately described with the inclusion of the cubic 
terms [Q. It was further shown that the induced cubic 
nonlinearity affects the phase modulation of cw waves, 
enabling efficient switching |17|, and that its strength 
can be increased by modulating the grating jist . 

For the more general QNPCs considered here, the in- 
duced cubic nonlinearity depends on both the DC part 
and the modulation part of the nonlinear grating, but is 
independent of the linear grating. The cubic terms may 
lead to either a focusing or a defocusing effect, depend- 
ing on the relative intensity of the fields and the sign 
of the phase mismatch /?, since sign(K) = sign(/3). The 
strength of the effective x''^'' nonlinearity depends on the 
difference in the linear grating strengths at the funda- 
mental and SH frequencies and on the DC component of 
the nonlinear grating. We thus recover the well-known 
effect that the interplay between the linear and nonlin- 
ear gratings can increase or decrease the effective x''^'' 
nonlinearity, depending on the physical situation 

The averaged model (^) has stationary, localised soli- 
ton solutions of the form wo{x, z)=e^^^w{x)/\p\ and 
vo{x, z)—e'^'^'^^v{x)/ p, which obey the equations 
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where 7=7/|/3p depends only on k, ai — 02, and da/d. 
The slowly varying approximation gives valid solutions 
when the soliton period is longer than the grating period, 
i.e. when the soliton parameter A is small, |A| <C \k\. 

It is important to stress that Eqs. (|^) cover a much 
more general situation than in [Q, which only consid- 
ered a nonlinear grating. A given value of the parameter 
7 represents a ra'nge of physical situatio'ns with different 
cornhi'natio'ns of li'near a'nd ■no'nli'near grati'ngs. In Fig. ^ 
we illustrate representative combinations for exact phase 
matching (/3=0) that all give the same value of 7. 



of Eqs. (0) using a standard relaxation technique. Bright 
soliton properties were investigated in [p^ , but not for 
normal dispersion, 7 > 0. In Fig. |^ we present the prop- 
erties for normal (7=0.02) and anomalous (7=— 0.02) dis- 
persion, together with the zcroth-order solution (7=0). 
The ratio R = v'^ (0) / 'w'^ (0) of peak intensities, shown in 
Fig. ^(a), confirms that the zeroth-order approximation 
becomes increasingly inaccurate for large A. Also, for a 
given 7, R approaches the same limiting value as A in- 
creases, regardless of the value of /3. In this limit the 
SH is stronger than the fundamental for 7 > {R > 1) 
and much weaker for 7 < (i? ~ 0). The total power 
P= X^(^^ + 'w'^)dx is shown in Fig. ^(b). For /3 > this 
reveals the interesting property that the power threshold 
for existence, decreases for 7 > and increases for 7 < 0, 
as compared to the zeroth-order value. 
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FIG. 3. QNPCs with the same value of 7. The linear (non- 
linear) grating is shown with a dashed (solid) curve. 

The first, simple case (/=!) is typical for domain in- 
version in ferroelectric materials, such as LiNbOs. It 
has only a nonlinear grating that flips from positive to 
negative with no DC component, do—0. The second 
case has a DC level of x'^) nonlinearity, do/d=3, but 
no linear grating, corresponding to the nonlinear part of 
the LiNbOs/HiLiNbOs structure reported in The 
third case is the GaAs/GaAlAs structure reported in , 
which has a nonlinear grating with do /d=5/3 and a linear 
grating with (oi— a2)//9=— 0.07. For the second and third 
cases to give the same 7 as the first, the grating wavenum- 
ber K must be multiphed by /=48.5 and /=26.6, respec- 
tively. Physically this can be done by changing the input 
beam width xq, thus maintaining (3=0. 
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FIG. 4. Soliton properties versus the internal parameter A 

for 7=0.02 (solid), 7=-0.02 (dashed), and 7=0 (dotted), and 

three values of the residual mismatch /3. (a) Ratio of peak 

intensities 7?={;^(0)/w^(0). (b) Total power P. 

We have numerically found the bright soliton solutions 
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FIG. 5. Ratio of peak intensities ii={i^(0)/i5^(0) versus A 
for K=10/ (continuous), k=4/ (dot-dashed), and k=— 10/ 
(dashed), and for the zeroth-order solution (dotted). The 
average of the propagating solitons are shown with /=1 
(crosses), /=48.5 (plusses), and /=26.6 (circles). /3=0. 



The bright soliton solutions of the average model (|5|) 
were tested for the three QNPCs of Fig. ^ by mapping 
them back to the variables w and and using them 
as initial conditions in simulations of the field Eqs. (|^). 
The evolution consists of small, regular oscillations su- 
perimposed on the slow average beam. Properties of the 
propagating solitons were calculated by averaging over 
an integer number of grating periods and were then com- 
pared with the predictions of the average model. Figure 
|5| displays the ratio of peak intensities versus A for exact 
phase matching, /3=0, and reveals that, for both anoma- 
lous and normal dispersion, the solutions of the average 
model are accurate for small A and large as expected. 
Even when |k|=4 the first-order solutions provide a much 
better fit than the zeroth-order solutions. Our analysis 
thus shows that bright solitons exist and propagate stably 
in QNPCs with many types and combinations of linear 
and nonlinear QPM gratings. 
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However, this is provided one is careful and does not 
eliminate the effective x^^^ nonlinearity by using a grat- 
ing with K—2{a2 — ai)do/d. This could happen in real- 
istic QNPCs without violating the assumption |k| 3> 1. 
For the LiNbOs and GaP/AIP structures given in 
at phase matching {k—P), the presence of the linear 
grating changes the effective x^^-* nonlinearity by a fac- 
tor of T—l + 2do{ai — a2)/(rf«;)=f .4 and 0.3, respec- 
tively. The linear grating thus adds constructively in the 
LiNbOa structure and destructively in GaP/AIP. In fact, 
modifying the nonlinear grating in the GaP/AIP struc- 



ture slightly to Xf?''=40pm/V (max.) and Xh^''=19pm/V 
(min.) eliminates the effective x^^' nonlinearity entirely. 

Conventional average models [ff-^ would predict that 
no soliton could exist with no nonlinearity, p=0. How- 
ever, in the model the induced cubic nonlinear- 
ity predicts that solitons should still exist as solutions 
of nonlinear Schrodinger equations. In the case when 
the SH is strong {vo/wq = \/5exp[i(/3 — X/2)z]), a 
family of bright solitons, wq—^JX/ (27) sech(V2Aa;)e*'^^, 
exists for normal dispersion (7, A>0), and a family 
of dark sohtons, wq = a/A/(47) tanh(Y^|A[a;)e*'^^, ex- 
ists for anomalous dispersion (7,A<0). With no SH, 
bright solitons, wq — ^2A/|7| sech('y2Ax)e*'^^, exist for 
anomalous dispersion (7<0,A>0), whereas dark solitons, 
^0 = v'l'^l/7tanh(-\/|A|a;)e''^^, exist for normal disper- 
sion (7>0,A<0). We test these solutions by mapping 
them back to the variables w and v and launching them 
as initial conditions in simulations of the field Eqs. (Ilj). 
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FIG. 6. (a) Bright and (b) dark solitons propagating in 
QNPCs with no effective x''^'' nonlinearity. Shown is the 
scaled intensity of the fundamental for |k| = 100, /3=0, and 
(a) 7 < 0, A=l and (b) 7 > 0, A=-l. The SH is zero. 

Figure ^ shows the evolution of a bright and dark soli- 
ton with no SH over a distance of 25 grating periods. The 
SH displays small, regular oscillations around the mean 
value zero, corresponding to the oscillations of the funda- 
mental seen in Fig. ||. The simulations thus confirm that 
bright and dark solitons can indeed propagate in a stable 
manner in QNPCs with no effective x*-^^ nonlinearity. 

In summary, we have shown that bright solitons exist 
and propagate in a stable manner in ID quadratic non- 
linear photonic crystals (QNPCs) with many types and 
combinations of linear and nonlinear QPM gratings. By 
deriving first-order averaged equations, we have shown 
that such QNPCs have an induced cubic nonlinearity, 



and we have numerically found previously unknown fam- 
ilies of bright solitons. Even with no effective quadratic 
nonlinearity, the QNPCs support both bright and dark 
solitons due to the induced cubic nonlinearity. We have 
found analytical expressions for these solitons and shown 
that they propagate in a stable manner. 

Dark solitons are always unstable in homogeneous x''^'' 
media in settings for type I SHG, due to modulational in- 
stability of the back-ground plane- waves . Our results 
show, for the first time, a dark soliton that appears to 
be stable under propagation, with the stabilizing mech- 
anism necessarily originating from the photonic crystal 
structure of the QNPC. An analysis of dark solitons in 
the general case is being carried out. 

The project is supported by the Danish Technical Re- 
search Council through Talent Grant No. 9800400. 
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